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Abstract —The Stirling heat engine operating in a closed regenerative thermodynamic cycle is 
analyzed. Polytropic processes are used for the power and displacement pistons. Following regener¬ 
ation, the maximum power density and efficiency are found and the compression ratio at maximum 
power density is determined. © 1997 Elsevier Science Ltd. 


INTRODUCTION 

The use of a regenerative heat exchanger (RHX) distinguishes the Stirling heat engine from the Carnot 
engine [1-7]. A power analysis of an endoreversible Stirling cycle was conducted previously and the 
maximum power (MP) and efficiency at maximum power (Tj* max ) based on higher and lower tempera¬ 
ture bounds were obtained for the cycle with perfect regeneration [8], The primary purpose of this 
study is to obtain a more realistic regenerator, which is obtained by introducing polytropic processes 
at the power and displacement pistons. By application of the MPD technique [9], maximum power 
density (MPD) conditions of the Stirling heat engine with realistic RHX are investigated. 


GENERAL THEORY AND MODEL 


A schematic of the equipment and the T-s diagram of the Stirling cycle for an ideal-gas working 
fluid are shown in Figs. 1 and 2, respectively. The exponents of the polytropic processes must be within 
the limits 0 < n, < 1 and 0 < n c < 1. The cycle compression and temperature ratio parameters are 
introduced as r v = v 4 /v 2 > 1 and t/r = T z /T 4 < 1, respectively. Using r v and i fi, the ratios of the tempera¬ 
tures at the ends of each cycle branch become 


A/T, = T 3 /T a = r*r\ TJT, = r'~\ A/A = r? 


(1) 


From Eq. (1), the net power of the Stirling cycle is 



Fig. 1. Schematic of equipment for the Stirling cycle. 
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Fig. 2. T-s diagram for the Stirling cycle. 

lV net = mRT 4 [n t ( 1 - - «,) - \- 1 )/( 1 - « c )]- 

After equating heats transferred at the regenerator, the regeneration condition 

i|r=(l -C "c- 1) 

is obtained. From Eq. (3), 

n c (ijj,n t ,r v ) = 1 - ln[l + (!/«/>)( 1 - r?r')]/ln r,„ 


Then, 


= mRT 4 ( 1 - rjt -1 )•!>, 


where 


( 2 ) 

(3) 

(4) 

(5) 


<t> = 


n t /(l - n,) - n c (4/,n,,r v )/[l - n c (i/r,n„r v )] 


( 6 ) 


It is found from Eq. (5) that n, > n c . From Eq. (3), r v < i//”, 1,1 is the requirement to ensure n, > n c . 
W d is obtained by dividing Eq. (5) by v 4 , viz., 

= W n Jv 4 = mRT, (r;' - nr 2 )<t>/v 2 . (7) 


The transferred heat is 


Qh = = m c v T 4 ( 1 - Cr' )(k - «,)/(1 - n t ). 


( 8 ) 


In view of Eq. (3), the thermal efficiency of the cycle obtained from Eqs. (5) and ( 8 ) is 

tj = (*-1)(1—«,)/(*-n t )<I>. (9) 


RESULTS 

Since and W dmix are known, we may compare ( W ncl /lV ma j and (W max / W dmax ) as functions of 
r). As may be seen in Fig. 3, the efficiency at the MPD ( 77 *. ) is greater than 17 * . 

d max max 
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Fig. 3. Comparison of work ratios vs. r\ (</>= 0.2, n, = 0.75). 


Figure 4 shows the effect of temperature on the two efficiencies tj*, and t} h>i , where we have 

max d max 

made use of t/». It is seen in Fig. 4 that the r\ w , values are higher than the ti* values at the same 

if/. Figure 5 shows that the r v ’s at the MP are always greater than those of the MPD. Stirling heat 
engines designed by considering the MP conditions will have larger dimensions than those designed 
at the MPD conditions. The difference between the sizes of Stirling engines designed by using these 



Fig. 4. Effects of cycle—temperature ratio ip on the thermal efficiencies at maximum power and maximum 

power density (n, = 0.75, k = 1.4). 
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Fig. 5. Variations of the cycle compression ratios at maximum power (vjv 2 )^ _ and at maximum power 

density (v 4 /v 2 )^ d with ijj (n t = 0.75). 


criteria is striking, especially for small <// values. It is observed from Fig. 6 that r\ increases with decreas¬ 
ing r v and increasing n t . 

The optimum working region of the cycle is shown in Fig. 7. The points (a) and (c) indicate the 
MP and MPD conditions, respectively. The condition of maximum efficiency occurs at point (b). There¬ 
fore, the arc (a)-(b) shows the optimum working conditions when power output is the main interest. 
On the normalized MPD curve, the arc (c)-(b) represents the optimum working conditions if the engine 
size and economic aspects are the main interests. The point (c ) has the advantage of having r\ corre¬ 
spond to T) u<i . The point (a) refers to Vw milK for the engine design in the MPD analysis. The point 

(e) is the intersection of the curves. When the working conditions at point (e) are determined, the 
engine has the following lower and upper power bounds: 


1 > (W n JW mM ) e > (W net .W max ) c , 
1 > OWc max )e > (WV< a > 

These two inequalities order the efficiencies as follows: 

> Ve > 


( 10 ) 

( 11 ) 


( 12 ) 


The arc (e)-(c') constitutes the optimum working region. The points (e) and (c ) show good agreement 
with the results of a recent study by De Vos [10] who stated that the optimal point for applications 
lies somewhere between the MP and Tj max points. Point (e) is important for high-power designs and 
point (c ) for engine size. 


CONCLUSION 

In the present study, we used polytropic processes in the power and displacement pistons to create 
a finite temperature difference at the regenerator for the purpose of obtaining a realistic RHX. For the 
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Fig. 6. The thermal efficiency rj vs n, for different cycle—volume ratios r,, (<p = 0.2, k = 1.4). 



Fig. 7. A detailed view of the optimal cycle region (i/i = 0.2, n, = 0.75). 
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determination of design parameters to reduce the engine size, the MPD technique was applied and the 
thermal design bounds of a Stirling heat engine determined. 
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NOMENCLATURE 


c v = Constant—volume specific heat 
k = Ratio of specific heats c p /c„ 
m = Mass-flow rate 
n c = Polytropic exponent at the 
displacement piston 
n, = Polytropic exponent at the power 
piston 

(2h = R ate °f polytropic heat transfer 
from the high-temperature heat 
source to the working fluid 
Q l = Rate of polytropic heat transfer 
from the working fluid to the low- 
temperature heat source 
R = Universal gas constant 
r,, = Compression ratio v 2 /v 4 
r 4 = Maximum cycle temperature at the 
power-piston exit 


T 2 = Minimum cycle temperature at the 
displacement-piston exit 
W d = Cycle power density W net /v 4 
= Maximum cycle—power density 
VV mlu = Maximum cycle—power output 
W„ cl = Net cycle—power output 
v = Specific volume 
v 4 = Maximum specific volume 
v 2 = Minimum specific volume 
ip = Temperature ratio T 2 IT 4 
tj = Thermal efficiency of the cycle 
rjmax = Maximum thermal efficiency 
T), d|nB = Cycle thermal efficiency at 

maximum cycle—power density 
tj^ = Cycle thermal efficiency at 

maximum cycle—power output 



